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Abstract 


The author's model for anisotropic solar cosmic ray propagation 
gives 2 coupled, partial differential equations for the intensity and 
anisotropy of solar cosmic rays propagating with finite speed V in 
an inhomogeneous medium. The model is used to study the effect of 
the solar shell on solar cosmic ray propagation. It predicts an 
exponential decay, regardless of the observer's position. It predicts 
that when the observer is near the center of the shell, to 

30, ( t D =decay time, t Q =onset time) and A m «157o, of t m /t Q ~3 to 5 
(t m =time of maximum), consistent with observations of relativistic 
particles on Feb. 23, 1S56. When the observer is between the shell 
and the sun, the model predicts that oscillations might be observed 
near maximum intensity. When the observer moves away from the sun 
and the shell, the propagation is diffusive, but there is an increasingly 
large presistent anisotropy which serves as a measure of the width 


of the shell. 
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I . Introduction 

In a companion paper, Burlaga (1969) presented a simple, 
phenomenological, 1-dimensional model for anisotropic solar cosmic 
ray propagation. This was used to study the general effects of a 
diffusing region near the sun. In the present paper the mathematical 
model is discussed further with emphasis on the limiting results for 
a continuous medium, and it is applied to examine qualitatively some 
basic effects of the diffusing region (solar shell) which is believed 
to extend from somewhere near the earth to a few or several AU. 
(Meyer, Parker and Simpson, 1956; Farker, 1963; McCracken et al. , 
1967; Fan et al., 1968; Lanzerotti, 1969; Burlaga, 1967). Whereas 
most of the earlier work emphasized the characteristics of the total 
cosmic ray flux (I) measured at 1 AU, che present work examines some 
of the qualitative results that might be seen by detectors which 
measure the anisotropy and intensity of solar particle fluxes over a 
wide range of distances from the sun. 
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II. The Model 

Basic Equations . The model is as described in detail by Burlaga 
(1969). We consider propagation in a 1-dimensional, semi-infinite 
medium with a source at x^=0. The medium is represented by an infinite 
number of point scattering centers which are equally spaced with 
separation \ and extend from the source to infinity. The basic 
equations are 


f L+l, T+l = P L f L,T f (1 ' P L )f L,T 


f, , m = (l-P T )f" + P T f; 
L-l , T+l 1/ L,T L L . i' 


( 1 ) 


where f + L T is the probability that a particle is moving away from 

the source after a collision at L-l, T-l and is just approaching 

point L after T-l collisions, f“ _ is the corresponding probability that 

a particle is moving toward the source, and P^ is the probability 

that a particle will be reflected when iu encounters a scattering 

center at x^. 

The Equations in the Limit of a Continuous Medium . Consider 
the limiting form of (1) when \- 0 . Since \ = Vt , where T is the 
time between successive collisions and V is the particle speed, the 
limit \-*0 implies that T-O. In this limit we can write ^l+ 1 T+l ~ 
fi +Vt , t T +T) which can be expanded in a Taylor's series to give 
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^L+l, T+l ~ f ^ X L ,t T^ + — ^ X L’ t P > V ' + **L* t T^ 

’ ox dt 

Putting this into (1), rearranging terms, and introducing the finite 

variable 


D(x) = TV 

= \V 

(2) 

2P(x) 

2P(x) 



gives 


-t- 

df 

+ v M + = 

_ v 2 

(f‘ - 

f ) 


dt 

dx 

2D(x) 




df" 

_ V of" = 

2 

V 

(f + - 

f") 

(3) 

dt 

dx 

2D(x) 





Note that in the limit T-*0, P(x) is on the order of r, since V and 
D(x) remain finite. Thus, in the process of taking the limit, D(x) 
has replaced P(x) as the basic function characterising the inhomogeneity 
of the medium. Adding and subtracting (3) gives the following equations 


31 dJ _ Q 
dt dx 




(4) 

J = - D(x) dl 
dx 

- D(x) dJ 
dt 


where I = f + +f 


(5) 

and J = (f + -f )V 


(6) 


Clearly, I is proportional to the total flux and J is proportional 
tc the net flux. From (4) one can show that 



dJ 

dt 




k [D(1 ° 


D(x) I? * 0 


dJ 

dt * 


dD(x) 

dx 


(7) 

( 8 ) 


In the limit V- 00 , (7) reduces to the familiar equation for diffusion 
in an inhomogeneous medium. In the case D = constant, (7) and (8) 
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reduce to the equations that Axford (1965) derived for anisotropic 
propagation in a homogeneous medium. When D = constant, (3) are 
identical to the equations which Fisk and Axford (1969) derived from 
the Boltzman equation for bi-directional scattering with P(x) = 0.5. 

The solutions of (7) and (8) with D = const, cannot give large anisotropies 
except at early times, as Fisk and Axford have noted. This is because 
there are scattering centers everywhere which always tend to reduce 
the anisotropy. 

Although it is customary to use solutions of the differential 
equations to test the theory, it is worth noting that the basic 
physics of the problem is contained in the expression for J in (4). 

It would be worth testing this directly, since it does not inv lve 
initial conditions, boundary conditions or assumptions about D(x) , 
and it gives a value of D(x) at the observation point which could be 
compared with the value computed from the power spectrum of the 
magnetic field. 

Numerical solutions of general equations for a homogeneous medium . 

Although P = constant is not very realistic for the interplanetary 
medium, it is worthwhile to consider solutions for this case because 
it is the simplest case and thus basic for an understanding of the 
model. 

When P = .5 the model reduces to the elementary random walk problem 
(Chandrasekhar, 1943). In this case one finds from numerical computations 
that 

I = 2 | exp(-L D 2 /2T) (9) 

where L Q is the number of scattering centers between the source and the 
observer and T is the total number of collisions that each particle has 
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experienced. This formula is valid only in the diffusion limit, 

L 0 «T and T large. In practice this is a remarkably good approximation, 

as Lord Rayleigh (1919) has noted, and valid shortly after particles 

first arrive at L . An improvement is obtained by using the solution 
o 

of the telegraph equation (Axford, 1965), but for a homogeneous 
medium the principle difference between this solution and (9) is that 
the telegraph equation gives 1=0 when T<L Q . 

When P<.5 we find that an equation like (9) describes the results 
obtained from the model in Section II, as long as P and L c are not 

too small. In particular, ijj is found that 

_ L o 

TFT 

1 = 2‘ 2 e , T> 1.2 L q (10) 

rrrT 


where 

r = (1-P)/F, 

and 

A s J/I = Is ± , T » T . (11) 

2P T 0 T 

Here T m is the number of steps corresponding to the maximum of I, and 
T 0 = L q is the number of steps between the source and the observer. 

Thus, apart from a constant multiplier, the intensity-time profile in 
a medium with P < .5 and a given X is equivalent to propagation in a 
medium with P=.5 and X' = rX. Axford (private communication) has 
pointed out that a similar result was obtained analytically by Goldstein 
(1951) for P=constant. 

To transform the above results to dimensional form, let 
x=LX , t-TT , and 


D = rXV 

2 


( 12 ) 
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Then 


and 


2 



D t m 1 
A = PrVx 7 


(13) 


(14) 


In the diffusion limit \-*0, v-°°, P-0, and Pr-*1; thus, with the 

2 

relation t = x_/(2D), one obtains 
m o ' ' * 

A = x o = .£ ££ (15) 

2Vt N dx 

Equations (13) and (15) are the usual results for diffusion in 


1 -dimens ion 
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III. Solar Shell 

The solar shell is the cosmic ray diffusion region in the 
vicinity of the earth. A thorough mathematical study of the shell 
requires a knowledge of P(x) , inclusion of the effects of the solar 
envelope, and consideration of both parallel and transverse diffusion 
(Burlaga, 1967, Bukata et al., 1969). Such an investigation involves 
several parameters which are not known experimentally, so a parameter 
study at this time is neither practical nor illuminating. Our approach 
will therefore be to study an idealized model which contains what we 
consider to be essential fe.itures of the real situation. The aim is 
to obtain a qualitative understanding of the effects of a shell which 
might be observed by spacecraft near and far removed from the earth's 
orbit. When more measurements are available, particularly measurements 
concerning P(x) (actually D(x)) which can be obtained by deep space 
probes, it will be meaningful to apply the model to fit actual 
observations in detail. 

We consider one-dimensional propagation as described in Section II. 
This will be a reasonable approximation if the flare is near a 
magnetic field line which passes through the observation point. Recent 
calculations of the diffusion constant based on magnetic field observations 
suggest that diffusion is essentially 1 -dimensional near 1 AU , (contrary 
to the assumption of Burlaga (1967)), but transverse diffusion may 
still occur in tha envelope and give rise to East-West effects in 
essentially the same way as calculated by Burlaga). It is well known 
that the effect of additional degrees of freedom in diffusive 
propagation is manifested by a power-law dependence in time. For 
example, the Green's function for 1-dimensional diffusion is 
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l 2 

proportional to t *exp(-x /4Dt) while for 3-dimensional diffusion it 

-3/2 ■? 

is proportional to t exp( -x~/4Dt) . Thus, the important qualitative 
features of the propagation will not be lost ay considering 1-dimensional 
propagation. 

As in previous work (Burlaga 1967, 1969), we shall assume an 

instantaneous point source and a non-absorbing boundary at the 

sun. We shall assume that the sun's surface reflects particles, 

but replacing the sun by interplanetary scattering centers would not 

qualitatively change the intensity-time profiles. To avoid specifying 

parameters characterizing the solar envelope, we shall neglect the 

envelope altogether. Sow, it is known that the envelope tends to 

store particles and thus acts as a modified source (Burlaga, 1969; 

Shishov, 1966). For the >1 bv particles on May 4, I960 the characteristic 

"injection time" for this source w-as t£^15 min. For the much lower 

energy particles on March 24, 1966, the same type of analysis gives an 

injection time t.j~fcO min. When the characteristic times of the observed 

intensity- time profiles are much greater than t,- , one can assume an 

instantaneous source with little error. When the characteristic 

times approach t^ the "smoothing" effects of the envelope will be 

seen, and when the characteristic times are >t the approximation is 

i 

invalid. 

The geometry of the solar shell is poorly understood. At times 

there seems to be little scattering between the solar envelope and 

the earth (McCracken et ai.,1967 ; McCracken, 1962) and in general the 

diffusion seems to cease beyond ~2 AU (Burlaga 1967). As a zerot'. 

2 

approximation to ?(x) we shall take Caxp(-(L-L { ,)/L v? ) . The assumption of 
symmetry is introduced solely for simplicity and has no basis in 
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observational facts. This is a point which must be examined by deep- 
space probes. The shell is centered at L g and contains ~2L W effective 
scattering centers. The characteristic diffusing length is determined 
by putting the observer at L Q . (Note that L Q is now independent of 
the solar envelope). For simplicity, we set C=.5. Equivalant 
qualitative results could be obtained for smaller C by increasing 
and L w , but the case C=.5 will serve to reveal the basic features of the 
propagation. 

With the above approximations, the model is rather simple. The 

results will n 'w be discussed for three different cases: L ~L , 

s o’ 

L < L , and L > L . 
so so 

L = L . First consider L = L , corresponding to an observer 
o s os ° 

in the middle of the shell. A typical intensity-time profile is shown 

in Figure 1, which gives the solution for L =20, L =20 and L =10. 

o s w 

The sunward flux rises gradually, reaches a maximum at T , and then 

decays exponentially with an e- folding time T^. The anti-sun flux does 

the same, with the same T^ but a larger T . The results can be summarized 

D m 

by the parameters T m /T o , T^/T o , A r( and A^,, where T^ is the onset 
"time" (the number of steps required to reach L^) , A m is the anisotropy 
at the time of maximum intensity of the total flux, and A m is the 
anisotropy at /ery large T. Here anisotropy is defined by the equation 


A 


f+ - f 
f + + f" ’ 


Several models were examined, with L Q =10, 20, 40, and ranging from 

3.33 to 40. It was found that the parameters describing the intensity- 


time profiles are determined primarily by L . Figure 2 shows 

w 


T / T , T_/T , and A 


m 


D o’ 


versus L for the various cases. Note 
w 


-1 
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the linear relationships. Also note that determines each of the 

3 observable quantities. Conversely, measuring these 3 quantities 

provides a strong test of the applicability of the model since they 

must all be related as required by the inferred L^. For example, 

T d /T^ = 15 implies that L„ = 10 which in turn implies that T m /T 0 =5.5, 

and A ^30%. 
m 

When the shell is very thin (1,^5) the profile is similar to that 

in Figure 1, except that there is a "precurser". This is caused by 

the relatively large number of particles which arrive at T Q without 

scattering because of the small L^. 

Observations of relativistic solar particles on Feb. 23, 1956, 

reported by Meyer et al. (1956), show an intensity-time profile which 

is similar to tnose in Figure 1 and a rapid approach to anisotropy 

as implied by the model for L 0 =L . Furthermore the flare was near 

the base of the field line passing through the earth (Burlaga, 1967), 

and the smallest characteristic time, t m ~130 min, was much larger 

than the time of escape from the envelope (15 min) for similar particles 

on May 4, 1960, so one might expect the above model to apply. Table l 

shows t m /t and t,./t for particles monitored at Ottawa (>1 bv) , 
m o Do ~ 

Chicago (>1.9 bv) and Wellington (>3.4 bv) . The values of t m and t 

are from Webber (1964) and the values of t are rectilinear propagation 

o 

times to e«rth for protons of the indicated rigidities. For the given 
t m /t one can find the L from Figure 2 and then a predicted t_/t , 

lu O y Do 

and A . Table 1 shows that the observed and predicted values of 
t /t Q are in reasonable agreement. (Exact agreement would be 
coincidental because of the assumptions in the model). The quantity 
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A was not measured directly, but the values shown in Table 1 
m 

are consistent with the observed rapid approach to near isot' ">py 
(Meyer et al. 1956, Lust and Simpson, 1957). Thus, the observations 
suggest that the middle of the solar shell was near the earth on 

Feb. 23, 1956. 

L < L . 
o s_ 

Now consider the case of an observer who is between the run and 

the center of the shell. Figure 3 shows f , f , and A versus T for 

L =20 and L =30. This illustrates the general case which shows a 
w s 

strong beam of particles arriving at T^ and moving toward the scattering 

region. Particles are reflected from the shell and mcve toward the 

sun were they are again reflected. The process continues so that one 

observes a damped train of oscillations in f and f which are out of 

phase such that the anisotropy also fluctuates. At late times the 

decay is exponential and the anisotropy tends to zero. The decay time 

is much larger than in the case L q =L s , and it depends appreciably on 

L as well as L . This is shown in Table 2. 
s w 

The 'period' of the oscillations is essentially constant. It 
equals 2(L -L ) when L <20 and is somewhat smaller for larger L . 

n s w w w 

The damping is simply the result of diffusion. Thus a thicker shell 

gives smaller amplitude oscillations and the amplitude of the oscillations 

is smaller when the observer is closer to the center of the shell. 

As one would expect, the oscillations are found to be very small when 

L 0 ~L w ; in chis case the shell extends to the sun and particles 

diffuse everywhere. 
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In order to observe these oscillations their period, which is 
approximately twice the time required for a particle to move from the 
sun to che shell, must be larger than the time for particles to escape 
from the envelope. If the shell center were at ~2 AU the period would 
be ~45 min. for 1 bv particles so one might expect to see the oscillations 
at the earth. Even if the shell were centered at ~1 AU, the period 
would be ~12 min. which is probably comparable to the ejection time 
for the envelope, and weak oscillations might be seen by instruments 
moving close to the sun. The real beam is not mono-energetic, as 
assumed in the model, and there will be complications arising from 
pitch angle scattering, mirroring, patches of turbulence, etc. All 
of these effects tend to damp the oscillations, so one cannot say 
from our model that the oscillations will be seen, only that they 
mi^ht be seen. If they are seen they can be used to probe the shell. 

Oscillations in solar cosmic ray fluxes have been observed (see 
the review by Williams, 1969) , but they are usually non-dispersive 
or have the 'wrong' period. (It is conceivable that the period could 
be velocity independent if the characteristics of the shell were 
appropriately velocity dependent, but this would require a special 
coincidence.) McCracken (1962) has reported particularly interesting 
observations of relativistic particles on Nov. 12, 1960 from several 
stations which view in various directions. Fluctuations in intensity, 
which could be interpreted as oscillations with a period ~20 min were 
observed at 2 stations which looked toward the sun, but the anti-sun 
flux d'^es not show oscillations with the proper phase. The very large 
anisotropies for the first 30 min. suggest that particles were stored 
for nearly this time in the envelope, which would also argue against 
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the detection of 20 min. oscillations. Thus, there is a s yet no 
convincing evidence for the oscillations predicted by the model. 

However, fhe shell is probably near the earth and, as noted above, this 
situation is not ideal for detecting the oscillations. They are more 
likely to be seen by instruments at~.5 AU. Relatively long decay 
times are likely to be observed at such distances even if the oscillations 
are not seen. 

Lp > L s . 

Finally, consider an observer who is beyond the center of the 

shell. Models with L^Lg+l^, where L o =10, 20, 40 and L w =L q /2 show 

that f 4 " and f rise gradually to a maximum and then decline exponentially. 

A typical profile is shown in Figure 4. Figure 5 shows that the 

ratio T /T varies linearly with L in essentially the same way as 

shown in Figure 2 for L Q =L W . However, the decay time is much larger 

when L =L + L_/2 than when L =L , and it depends on L_ as well as 
o s o os 3 

L . The anisotropy is large, ranging from 307, when 1^=20 to 807„ 

when L =3.3, and A^A^ Figure 5 shows that A varies linearly with 

L . When L n =L_ + aL where a> 1, the anisotropy increases because 
w 3 w 

there is less backscattering beyond the observer. 

The observations at 1 AU never show a large anisotropy at late 
times, although small anisotropies have been observed (McCracken et 
al. 1967), and attributed to effects other than geometry. However, when 
instruments are carried toward the outer planets, large anisotropies 
at late times and large decay times should be observed. Clearly, 
measurements of the anisotropy would provide a means of determining 
the relative position of the shell and are thus highly recommended. 
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TABLE 1 . Parameters for the February 23, 1956, Event 



t m^ t o 

L 

w 

V'o 

L D y/t o^ t * 1 '^ 

A m (th.) 

1 bv 

5.2 

10.5 

21 

27 

-15% 

1.9 bv 

4.4 

9 . 3 

21 

24 

15% 

3.4 bv 

3.8 

8.5 

15 

22 

15% 
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TABLE 2 . Decay 'time' versus L and L_ for L < L 

w s so 


L 

L 

T 

w 

s 

D 

3.33 

13.33 

130 


16.33 

165 


20 

200 

5 

20 

210 


30 

285 


40 

260 

10 

30 

805 


40 

1095 


50 

1390 

20 

40 

2020 


60 

3140 


80 

4280 


100 

4635 

40 

80 

5680 


120 

7785 


160 


9200 
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Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


F IGURE CAPTIONS 

'Intensity-time' profiles for 'sunward' (+) and 
'anti-sun' (-) 'fluxes' for an observer at the center of a 
shell who?* thickness is comparable to the distance 
between its center and the source. Specifically, 
f(L,T) is the probably of finding a particle at the 
scattering center from the source after it has 
made T collisions. Note the gradual rise to maximum, 
the exponential decay, and the small persistent 
anisotropy at late times. The anisotropy is shown 
at the bottom of the figure. 

For events with L 0 =L w , such as that in Figure 1, 

T /T , T_/T and A are determined primarily by L 
m o D ° m w 

and vary with L w as shown here. The subscripts m, 

D and zero refer to the times of maximum, decay, and 
onset, respectively. 

When the observer is between the shell and the source, 
the model predicts oscillations in the fluxes and 
anisotropy, as illustrated here. Note that the shell 
has the same thickness and the observer is at the 
same position as the case in Figure 1 , but the shell 
is farther removed from the source. 

When the observer is beyond the shell, the 'intensity- 
time' profiles are 'diffusive' , but a large anisotropy 
is predicted dun g the decay because there is little 
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Figure 5. 


backscattering beyond the shell. 

When the observer is beyond the shell, T /T and 

in ° 

the anisotropy at maximum intensity and during the decay 
are linearly related to L w as shown here. 
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